function [idiosyncratic_draw] = draw_idiosyncratic_Simple(idiosyncratic,q,rho,...
    Sigma_v,shock)

for i = 1:size(idiosyncratic,2)
    % Find the position of the last non NaN
    
    start = find(~isnan(idiosyncratic(:,i)), 1, 'last')+1;
    finish = size(idiosyncratic,1);
    
    for t = start:finish
    idiosyncratic(t,i) = 0;
        
        for j = 1:2
            idiosyncratic(t,i) = idiosyncratic(t,i) + idiosyncratic(t-j,i)*rho(i,j);
        end
        
        if shock == 1       % Either compute the point estimate or shock it
        idiosyncratic(t,i) = idiosyncratic(t,i) + mvnrnd(0,Sigma_v(i));
        end
    end
end



idiosyncratic_draw = idiosyncratic;